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ABSTRACT 

There has been much debate about the origin of the diffuse 7 -ray background 
in the MeV range. At lower energies, AGNs and Seyfert galaxies can explain 
the background, but not above ~0.3 MeV. Beyond ~10 MeV blazars appear 
to account for the flux observed. That leaves an unexplained gap for which 
different candidates have been proposed, including annihilations of WIMPS. 
One candidate are Type la supernovae (SNe la). Early studies concluded that 
they were able to account for the 7 -ray background in the gap, while later work 
attributed a signihcantly lower contribution to them. 

All those estimates were based on SN la explosion models which did not 
reflect the full 3D hydrodynamics of SNe la explosions. In addition, new 
measurements obtained since 2010 have provided new, direct estimates of high-z 
SNe la rates beyond 2 ; ~2. We take into account these new advances to see the 
predicted contribution to the gamma-ray background. 

We use here a wide variety of explosion models and a plethora of new 
measurements of SNe la rates. SNe la still fall short of the observed background. 
Only for a £t, which would imply ~150% systematic error in detecting SNe la 
events, do the theoretical predictions approach the observed fluxes. This fit is, 
however, at odds at the highest redshifts with recent SN la rates estimates. 
Other astrophysical sources such as FSRQs do match the observed flux levels in 
the MeV regime, while SNe la make up to 30-50% of the observed flux. 


Subject headings: gamma rays: difuse background; supernovae; galaxies: active 
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1. Introduction 

The cosmic gamma-ray background is diffuse and its origin is diverse and remains 
partly unknown at various energy ranges. On the low energy side, from X-ray energies up 
to around 0.3 MeV, AGNs and Seyfert galaxies provide most of the emission (Madau et 
al. 1994; Ueda et al. 2003). At if ^ 0.3 MeV, the spectrum of AGNs and Seyfert galaxies 
sharply cuts off. From 50 MeV to the GeV range, blazars seem to be responsible for the 
observed flux (Zdziarski 1996; Sreekumar et al. 1998). However, the latest results from 
Fermi, in the GeV range, which show a higher gamma-ray background at GeV energies than 
previous results from EGRET, have called into question the former attribution to blazars as 
the main source (see discussions in Lacki, Horiuchi & Beacom 2014). Nevertheless, recent 
examinations of this issue hnd that at energies > 100 MeV, blazars account for ~50% of the 
background, while the other half is contributed by star-forming galaxies and radio galaxies 
(Ajello et al. 2015; Di Mauro & Donato 2015). 

The measurements in the MeV range have been provided by various space missions. 
The hrst exploration of the region between 1 and 5 MeV was made by the APOLLO 15/16 
missions (Trombka et al. 1977). The reanalysis of the Apollo data, the measurements from 
HEAO-A4 (Kinzer et al. 1997), the Solar Maximum Mission (Watanabe et al. 1999a), and 
GOMPTEL (Kappadath et al. 1996; Weidenspointner 2000) provided the basic empirical 
results on the diffuse gamma-ray background, in the range from 100 keV to 10 MeV. The 
slope of the emission spectrum exhibits a steep decrease with increasing energy, from a few 
hundred keV to 10 MeV, changing to a flatter slope around 10 MeV and beyond, revealing 
the need of an intense extragalactic sonrce in the MeV window. In this range of energies, 
the discussion of the origin of the backgronnd was revived at the beginning of the XXIst 
century and continued up to the present time (see Lacki, Horiuchi & Beacom 2014). 


Much of the current discussion centers on the possibility that Type la snpernovae are 
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able to produce the observed flux in the 0.3—3 MeV range, Ailing the gap between Seyferts 
and blazars. This possibility was first suggested by Clayton and Silk (1969), and further 
studied by The et ah (1993) and Watanabe et al. (1999b) (hereafter W99b). At those 
times, there were no empirical rates of SNe la available, in particular the function Ria{z), 
was not measured at the required high redshifts. Therefore, the authors had to estimate 
fluxes on the basis of the better known star formation rates in galaxies and a convolution 
with an assumed delay time distribution between the birth of the progenitor binary systems 
and the explosions. Depending on the delay time since birth of the system, the SNe la rate 
could be either too low or just right, to account for the measured level of the gamma-ray 
background (W99b). Ruiz-Lapuente, Casse & Vangioni-Flam (2001) (hereafter RCVOl) 
performed a study of the SNe la contribution in the MeV region, by considering a wide 
range of star formation rates in galaxies and different efficiencies in the production 
of SNe la, for each star formation rate. Ruiz-Lapuente (2001) found that not much under 
1000 Mq going into star formation gives rise to one SN la (720 ± 250 Mq produce 1 SNIa). 
The efficiency, defined as 


StJz) = 


Riaiz) yr ^Mpc 


.-3 


( 1 ) 


P:t.{z) Mq yr~^Mpc~^ 

was 1.41 ± 0.35 10 Mq“^ hgg. In that work, they found that the delay time between 
star formation and supernova explosion could shift the estimate up or down, but that most 
star formation histories would average out the effect, since the distribution of SNe la delay 
times appeared to be broad. The result of this study was that SNe la yield a background 
emission, in the MeV range, that can explain the extragalactic emission measured by 
COMPTEL and SMM. 


Just shortly afterwards, the topic was addressed by other authors (Iwabuchi & 
Kumagai 2001; Ahn, Komatsu & Hoflich 2005 —hereafter AKH05—; Strigari et al. 2005), 
and was more recently considered again by Horiuchi & Beacom (2010, hereafter HB2010). 
The latter authors obtained a range of possibilities concerning the SN la contribution, but 
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found that they are in general not able to account for the MeV background. The earlier 
’’negative” results can be traced to the hrst high- 2 : SNe la rates obtained by Dahlen et ah 
(2004, hereafter D04) which peak at a smaller 2 than current rates (Okumura et ah 2014; 
Rodney et ah 2014; Graur et ah 2014). However, HB2010 do not use the Dahlen et ah 
(2004) rates, but a common £t to data availablle up to 2010. These rates will be compared 
with most recent measurements. 

Both W99b and RCVOl (see also similar findings by The et al. 1993 and by Zdziarski 
1996) had concluded that there was room for SNe la to account for the unexplained MeV 
background. They did so by integrating the full gamma-ray spectrum emerging from SN 
la models during the early period of significant production by radioactive decay. In this 
context one includes the production of gamma-ray photons as well as their transport in 
the expanding ejecta. AKH05 take a SNIa rate (from D04) that is one order of magnitnde 
smaller than the one adopted by W99b. This may explain most of the discrepancy in the 
results. HB2010 calculations take the continnnm plus line 7 -ray spectra from the W7fm 
model (fnlly mixed W7 model) shown in W99b and the 5p0z22.23 model shown in AKH05 
(2005) (private commnnication). We find agreement with their results when using the same 
SNe la rates as HB2010 did. 

One interesting featnre of the SNe la contribntion, which is shared by the calculations 
presented in W99b, RCVOl, AKH05, HB2010 is that the predicted SNe la gamma-ray 
backgronnd contribntion has a spectral shape (power law in energy) that runs parallel to 
the observed fluxes up to 3 MeV. This match in slope is one of the primary reasons for 
considering SN la as a significant contribntor, and to focus on the disagreement with the 
flux normalization. 

Dne to the recent appreciation of SN la as cosmological probes, a significant increase 
in dedicated observational programs prodnced a greatly improved knowledge of their rates. 
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for a wide range of redshifts (some of the most recent contributions were mentioned above). 
We also know well the cosmological parameters required for the calculation {Hq, Qm, ^a)- 
Moreover, for the hrst time, a nearby SN la (SN 2014J) has been detected in the MeV range 
(with INTEGRAL: Churazov et al. 2014; Diehl et ah 2014). The latter authors account for 
the gamma-ray data with a white dwarf explosion having a small amount of ^®Ni at the 
outskirts and around 0.5 Mq in the innermost core. This model is similar to those used in 
the previous calculations of the gamma-ray background (and also to those that will be used 
here). Thus, we now have a hrmed up grasp of the gamma-ray flux emitted by SNe la. 

These developments led us to once more calculate the gamma-ray background from 
SNe la, now based on improved knowledge of the theoretical ingredients. This paper is 
organized as follows. In Section 2, we discuss the different numerical methods used to 
compute the gamma-ray fluxes from given SNe la models. The models used for the present 
study are described in Section 3. In Section 4, we present an update of the SNe la rates 
at various ^ and show the interpolations used for the background calculations. In Section 
5, we explain the gamma-ray background formulation and present the results. Section 6 
discusses the results and the hnal Section contains our conclusions. 


2. Gamma—ray escape 

In order to derive the emergent spectrum predicted by supernova models, we use the 
Monte Carlo code described in Burrows & The (1990) and The, Burrows & Bussard (1990), 
modihed to include Bremsstrahlung X-ray production and the iron fluorescence line at 
~6.4 keV (Clayton & The 1991; The, Bridgman & Clayton 1994). This code follows the 
gamma-ray emission of radioactive nuclei such as "^^Ti, ^®Ni, ^®Co, ®^Co. The code was used 
in the predictions by W99b. 
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In RCV2001, a modified version of the Ambwani & Sutherland (1988) code was used, 
which was also employed by Lehoucq, Casse & Cesarsky (1989) to study the radioactive 
output of SN1987A. It was further modihed (Ruiz-Lapuente et ah 1993) to include 
positronium formation in the ^®Co decay, giving rise to two-photon (parapositronium) 
and three-photon decay (orthopositronium), the energy distribution being that derived 
by Bussard, Ramaty & Drachman (1979). The Monte Carlo routine was then used for 
the study of gamma-rays from SNe la (Ruiz-Lapuente et ah 1993). A comparison of the 
predicted gamma-ray emission, for model W7 (Nomoto, Thielemann & Yokoi 1984), with 
that from W99b shows the gamma-ray fluxes predicted from the code used in RCV2001 to 
be 10% higher in the MeV regime than those from W99b. Differences of the gamma-ray 
yields resulting from different gamma-ray transport codes were studied in detail by Milne 
et ah (2004), who show that variations of ~10% can arise from the different physics and 
methods employed. With the use of the same code as in W99b, Burrows & The (1990), and 
The, Burrows & Bussard (1990), we are on the conservative side concerning the predicted 
gamma-ray emission (see Milne et al. 2004). 


3. Input models 

A wider range of predicted emerging fluxes in the MeV domain is obtained when using 
different input models for the SNe la, as recently shown by The & Burrows (2014). Looking 
at these results, one notices that there are ~20% variations, in the MeV range, between 
various input models. Here we pick hve different models in an attempt to delineate the 
upper limit constraint on the gamma-ray emision from normal SNe la. These models are 
W7 (Nomoto, Thielemann & Yokoi 1984), the fully mixed W7 model (W7fm) used in W99b 
(both models synthesize 0.58 M© of ®®Ni), the W7dt model of Yamaoka et al. (1992) (this 
one predicts a larger escaping gamma-ray flux than most of the other models explored by 
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The & Burrows 2014 as it synthesizes 0.78 Mq of ^®Ni, a quantity slightly larger than in a 
normal SNla), the spherically-averaged 3D model NlOO of a delayed detonation (Ropke et 
al. 2012; Seitenzahl et ah 2013), and the (also averaged) 3D model of the violent merging of 
two WDs of Pakmor et ah (2012). Those two 3D averaged models synthesize, respectively, 
0.604 Mq of ®®Ni in the model NlOO and 0.62 Mq in the averaged 3D violent merging of 
two WDs, quantities within the canonical range of 0.55-0.65 Mq of ^®Ni for normal SNe la. 

The first three models have been widely discussed in the literature. Their gamma-ray 
emissions are compared in Figure 1. The 3D model NlOO of Ropke et al. (2012) corresponds 
to an initially isothermal WD, made of C and O in equal parts, with a central density of 
2.9 X 10® g cm“^. It was ignited in 100 sparks (hence its name), placed randomly in a 
Gaussian distribution within a radius of 150 km from the WD’s center. After an initial 
deflagration phase, a detonation was triggered at every location of the flame where the 
criterion for a deflagration-to-detonation transition of Ciaraldi-Schoolmann et al. (2013) 
was met. 

The model from Pakmor et al. (2012), also a 3D model, simulates the violent merger 
of a 1.1 Mq and a 0.9 Mq WD. As the material of the tidally disrupted secondary WD hits 
the primary WD, a hot spot forms which leads to the ignition of C burning. The conditions 
reached there are sufficient to trigger a detonation that burns the merged object and leads 
to the explosion. 

The spectral evolution of the gamma-ray emission from the two 3D models above is 
shown in Fig. 1 of Summa et al. (2013), together with the spread due to different viewing 
angles for the maximum-light epochs in gamma rays of both models. Here, however, we use 
the spherically-averaged models, with the same treatment of gamma-ray escape as for the 
three ID models (W7, W7fm, and W7dt). Their comparison with W7fm is shown in Figure 


2 . 
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Ruiter et al. (2013) have shown that sub-Chandrasekhar pure detonation models 
can reproduce the observed peak-magnitude distribution of SN la (Li et al. 2011). The 
brightness of the explosion is mainly determined by the ^®Ni mass synthesized in the 
detonation of the primary WD. Sim et al. (2010) derived a relationship between the 
mass of that WD and the expected peak bolometric brightness. Ruiter et al. (2013) 
use the population synthesis data from Ruiter et al. (2011) to derive the theoretical 
peak-magnitude distribution. 

Here we estimate the contribution of this channel to the gamma-ray background, based 
on the fact that the distributions of WD and ^®Ni masses peak at 1.1 Mq and 0.6 Mq, 
respectively. In Figure 3 we compare the gamma-ray spectrum of a representative model 
of such mergings with that of model W7fm. The major difference between models is in the 
X-ray domain between 1 and 50 keV. However, in the range of MeV, the difference between 
these models that reproduce the bulk of SNe la in the optical is of order 5%. Fig.4 shows 
the difference between the models in different energy ranges. 

What has been shown in this Section is that the background arising from different 
explosion models, associated with different progenitor origin channels (but all aiming at 
explaining normal SNe la) induces only a small variance. The convolution of the gamma-ray 
yield resulting from several models, here exemplihed by model W7fm, with the supernova 
rates is explored in Section 5. 


4. SNe la rates 

Several groups have obtained SNe la event rates using either ground-based facilities, 
the combination of ground-based and space facility such as the Subaru/XMM-Newton 
Deep Survey (Okumura et al. 2014 being the most recently published search of this kind). 
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Fig. 1.— The total X-ray continuum and gamma-ray line fluxes per ®®Ni nuclei, integrated 
over the first 600 days. The green line shows the corresponding fluxes from the calculations 
of The et ah (1993) for the W7 model (Nomoto, Thielemann, & Yokoi 1984), the blue line 
shows the fluxes for the fully mixed W7 model (W7fm) used in W99b, and the red line shows 
the fluxes for the W7dt model from Yamaoka et al (1992), (see also The & Burrows 2014). 
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Fig. 2.— The total X-ray continuum and gamma-ray line fluxes per ®®Ni nuclei, integrated 
over the first 600 days comparing the calculations for the fully mixed W7 model (W7fm) 
used in W99b (blue line), for the (spherically averaged) delayed-detonation 3D model NlOO 
of Ropke et ah (2012) (green line), and for the (also averaged) 3D model of the violent 
merging of two WDs of Pakmor et ah (2012) (red line). 
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Fig. 3.— Total time-integrated X-ray continuum and gamma-ray line fluxes per ^®Ni nucleus 
, comparing the calculations for an average model of the sub-Chandrasekhar pure detonation 
models (Sim et ah 2010; Ruiter et ah 2013) (red line), and for the W7fm model of W99b 
(blue line). 
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Fig. 4.— The ratio of number of emergent photons relative to model W7. The bar graph 
shows the ratio of integrated number of emergent photons of models considered here relative 
to W7. Note that photons in energy bands below 0.6 MeV are not relevant to the CGB near 
MeV region, which is our main consideration, but we plot them for completeness. Please 
note the good agreement within 5 % of the 3D models with the w7fm model as well as with 
the W7 model. The model w7dt has more ^^Ni mass than a normal SNIa. 
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or the Hubble Space Telescope (Rodney et ah 2014, and Granr et al. 2014 have presented 
the most recent resnlts from snch searches). One goal of those surveys is to gauge the 
delay-time distribution (DTD) of SNe la, i.e., the time that the SNe la progenitors spend 
between formation and explosion. 

Figure 5 shows how the different hts to compilations by Okumura et al. (2014), Granr 
et al. (2014) and Rodney et al. (2014) pass through the data up to z ~ 1 with similar lines, 
but they differ signihcantly at larger redshifts. For a close look at low and high z, we plot 
in Figure 6 all measurements made up to now, taken from the compilation in Table 4 of 
Granr et al. (2014) (superseded ones excluded)t 

The estimates of the SNe la rates are derived in different ways. Okumura et al. (2014), 

®In order of increasing redshuft, they are from: Gappellaro et al. (1999), Li et al. (2011a), 
Dilday et al. (2010), Dilday et al. (2008), Dilday et al. (2010), Granr & Maoz (2013), Blanc 
et al. (2004), Hardin et al. (2000), Dilday et al. (2010), Rodney & Tonry (2010), Perrett 

et al. (2012), Horesh et al. (2008), Dilday et al. (2010), Dilday et al. (2010), Perrett et 

al. (2012), Botticella et al. (2008), Dilday et al. (2010), Rodney & Tonry (2010), Perrett 
et al. (2012), Granr et al. (2014), Rodney & Tonry (2010), Perrett et al. (2012), Granr et 
al. (2014), Rodney & Tonry (2010), Perrett et al. (2012), Tonry et al. (2003), Dahlen et 
al. (2008), Pain et al. (2002), Rodney & Tonry (2010), Perrett et al. (2012), Melinder et 
al. (2012), Rodney & Tonry (2010), Perrett et al. (2012), Granr et al. (2011), Rodney & 

Tonry (2010), Perrett et al. (2012), Barbary et al. (2012), Dahlen et al. (2008), Rodney & 

Tonry (2010), Perrett et al. (2012), Granr et al. (2014), Rodney & Tonry (2010), Perrett 
et al. (2012), Rodney & Tonry (2010), Perrett et al. (2012), Barbary et al. (2012), Dahlen 
et al. (2008), Granr et al. (2011), Barbary et al. (2012), Granr et al. (2014), Dahlen et al. 
(2008), Granr et al. (2014), Dahlen et al. (2008), and Granr et al. (2011), and are shown as 
open circles in the Figure. 
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following previous work by the Supernova Cosmology Project, select objects having a flux 
larger than 5 at, where at is the background fluctuation within an aperture of 2 arcseconds 
in diameter. These objects are classihed as transients. The authors require that objects 
show at least a 5crf, increase in 2 or more epochs, in the i'-band. This requirement is 
accounted for in the rate calculation when computing the control time. For the supernova 
hosts for which spectroscopic redshifts are not available, they use photometric redshifts 
from the probability distribution function of the host galaxy redshifts. They discriminate 
against AGN using X-ray data from observations with XMM-Newton. Template light 
curves help to discriminate SNe la from SNe Ib/II. 

The hnal count of SNe la detected leads to the number of supernovae per comoving 
volume unit, rv{z). The number of SNe la expected in a redshift bin {zi < z < Z 2 ) is given 
by (Okumura et al. 2014): 

CTiz) 0 

Nexp{zi < Z < Z2) = / rviz)— - —V{z)dz (2) 

j z\ L z 47r 

where V{z)dz is the comoving volume in a redshift slice of thickness dz, 0 is the solid angle 
observed in the survey (in units of steradians), and CT{z) is the observer frame “control 
time”, i.e. the total time for which the survey is sensitive to a SN la at redshift ; 2 . The SNe 
la rates per comoving volume of the Subaru/XMM-Newton Deep Survey (Okumura et al. 
2014) are represented by blue triangles in Figure 5 and 6 together with the compilation of 
rates at all different z. They include a generous 50 % error in the systematic uncertainties 
due to dust obscuration leading to lack of detection of SNe la. In the hgures, we also 
plot the la upper and lower bounds of Okumura et al.(2014) rates as two blue long-short 
dashed lines. They represent the total, statistical and systematic, uncertainties of the points 
measured by Okumura et al. (2014). 

Graur et al. (2014) presents a compilation of all previous measurements of SNe la 
rates and new rates obtained by the GLASH collaboration in the range 1.8 < z< 2.4. The 
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procedure they use to compute the rates has aspects in common with Okumura et ah (2014) 
approach, as they look for candidates in this case with a flux 3a above the background and 
they rely on a good knowledge of the redshift of either the SN or its host galaxy. But it also 
includes some Bayesian approach in the classihcation of the SN with a probability obtained 
by comparing the light curve of the supernova with light curve models, as in Rodney et ah 
(2014). 

Rodney et ah (2014) (for the CANDELS collaboration) use technique and requirements 
different from those of Okumura et al. (2014). They use a Bayesian photometric 
classihcation algorithm based on the light curves of high-redshift supernovae. The main 
uncertainties in this approach are: the redshift-dependent prior describing the relative 
fraction of SN that are SNe la, and the assumed distribution of dust extinction values. 
The host-based prior comes from the algorithm galsnid of Foley and Mandel (2013), 
which exploits the relationships between SN and their host environments. They adopt 
a distribution of host galaxy extinctions for the dust extinction that then enters as a 
further probability into the SNe la classihcation algorithm. The mid-rate plus mid-dust 
combination gives their baseline classihcation probability, which determines how much each 
individual SN contributes to the total count of observed SNe la. The upper and lower 
bounds make the error bars. The SNe la rates from the CANDELS collaboration are 
plotted as green squares in Figures 5 and 6. 

The utilization of informed Bayesian priors to estimate the SNe la rate as a function of 
redshift was published earlier by Graur et al. (2001). Kuznetsova et al. (2007) had followed 
a Bayesian approach as well, to determine the rates in four redshift bins. 

The most conservative SNe la rate error bars derived up to date are those from 
Barbary et al.(2012). We show them as inverted black triangles in Figure 5 and 6. The 
study considers diherent dust models in estimating SNe la rates where it gives the highest 
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estimates of the observed SNe la rates, and also the largest uncertainties. The large error 
bars also result from the large uncertainties of dust extinction estimates. 

The D04 data are shown as red circles in Figure 5. The £t to these data appears as 
the red long-dashed line in the Figure. As it can be seen, we have fitted the D04 data in 
a similar way as in AKH05. We have considered, as in that work, that the first errorbar is 
reduced by taking the mean value of the upper and lower error bars of the hrst 2 :. 

One controversial aspect in the Icr and 2cr error bars from Okumura et ah (2014) is the 
way they are calculated. If we were to take the Icr upper and lower limits to the SN la rates 
by propagating the errors in the fitting parameters given by their expressions (10)-(12), 
the upper Icr limit would be at odds with the data from all other surveys, and also well 
beyond that indicated by the uncertainties of their own data. There is no questionable 
deduction, however, of the Icr errors of their individual measurements presented for first 
time in their 2014 paper. Nevertheless, the extrapolation of the htting formula beyond their 
data departs significantly from what has been found by other authors (Rodney et al. 2014; 
Graur et al. 2014). We keep their fitting formula, given that discrepancies start beyond the 
redshift range from within which most of the background arises, and we use the Icr and 2cr 
uncertainties as derived from their data points. 

To date there remain divergences in the measured SN la rates, despite the considerable 
increase of the data base in recent years. They reflect the differences in the procedures 
followed to derive the rates from the observations. In particular, how the SNe la rates 
estimate accounts for dust obscuration. 

The approach by Okumura et al. (2014) is well tested. It clearly discards SNe la from 
AGNs and SNe of other types. Gontamination by SNe Ibc is estimated and taken into 
account in the systematic uncertainty. Thes authors acknowldege the significant uncertainty 
related to dust extinction (50 %). Note the present disparity in the treatment of dust 
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obscuration and the very different approaches taken by different authors. 

We just report the published rates and their corresponding gamma-ray distributions 
based on three different estimates of these rates at the highest available z to show the 
level of uncertainty at high z while they are consistent at low z, where these three major 
estimates take into account all previously measured SNe la rates (see Figures 5 and 6). We 
calculate the gamma-ray background contribution corresponding to those SNe la rates. 
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Fig. 5.— SN la rates measured at different redshifts and diverse fits to them. Red 
circles are D04 data from Dahlen et al. (2004); green squares come from the CANDELS 
program (Rodney et al. 2014); blue triangles are from the Subaru/XMM-Newton Deep 
Survey (Okumura et al. 2014); inverted black triangles correspond to the HST Cluster 
Supernova Survey (Barbary et al. 2012). The red long-dashed line is a fit to the D04 data, 
done as in AKH05. The green, short-dashed line, is the same fit to the CANDELS data 
shown in Fig. 11 of Rodney et al. (2014), while the magenta dot-dashed line corresponds to 
the rates used by HB2010 (private communication). The solid black line is a £t to the rates 
given in Table 6 of Graur et al. (2014) (see next Figure for details), and the two black dotted 
lines are its Tier limits. Finally, the blue solid line is the same fit as in Fig. 11 of Okumura 
et al. (2014) to the Subaru/XMM-Newton Deep Survey data, and the two blue long-short 
dashed lines its Tier limits. Note that the upper la limit just fits the two highest measured 
rates. The ila regions around the rates from Okumura et al. (2014) have been traced by 
calculating the error using the statistical and systematic errors of the points measured by 
these authors. Note that the upper systematic errors include 50% extinction by dust, which 
explains why the upper error line is more distant from the main fit than the lower one. 
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5. Background contributions 

Based on the gamma-ray spectra in Figures 1-3 (number of escaping gamma photons 
per keV and per nucleus of ^®Ni synthesized in the explosion, as a function of energy E), 
calculated for different SN la models, we compute the average luminosity 1^{E) (photons 
s“^ keV“^) of a SN la from the amount of ^®Ni in the model and the time along which 
the SN has been emitting. In the input spectra, all the photons emitted from the time 
of the explosion to 600 days later were collected. Therefore, the number of active SNe la 
per unit of comoving volume, at any time, is that of those produced during the preceding 
time interval: R'ja{z) = const, x Rja{z), the latter being the comoving SN la rate (SN yr“^ 
Mpc“^). Thus, the contribution to the gamma-ray background of the shell at comoving 
radius r and with thickness dr would be: 

dL^{E, z) = ATiR'j^{z)l^{E)dV{z) (3) 

where 

dV{z) = dl^{z)d{dM) (4) 

dM being the proper motion distance (in Mpc). The flux received from that shell (in 
photons cm“^ s“^ keV“^) will be: 

dE^{E, z) = + 1)E, z] (5) 

di being the luminosity distance (in cm). The factor [z + 1), multiplying E, accounts for 
the redshift of the photons. We thus have: 

dF-,(E. z) = R',..(z)l,l(z + l)E|^£Md(d„) (6) 

Due to time dilation, there should be a factor (1 + z)~^ multiplying the comoving SN rate, 
but it is cancelled by the {1 + z) factor accounting for compression of the energy bins. 
Then, since di = {I + z)dM- 

= (X I + l)E]d{dM) 


dF.y{E, z) 


( 7 ) 
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Fig. 6.— Different fits to the currently available SN la rates data, up to a redshift z ~ 2.25. 
The solid black line is a fit to the rates given in Table 6 by Graur et ah (2014) and the two 
black dashed lines its ±la limits. The blue solid line running up to the upper edge of the 
diagram is the fit by Okumura et al. (2014) to the Subaru/XMM-Newton Deep Survey data, 
and the two dotted lines the +-la limits. The green solid line is the £t to the CANDELS 
data (shown as green squares) alone, from Rodney et al. (2014), while the magenta solid 
line corresponds to the SN la rates adopted by HB2010 (private communication). The open 
circles are the data from different surveys, covering from z = 0 upto2:~1.7 (see text for 
the references). Note, when comparing with the previous Figure, the vertical scale is here 
linear instead of semilogarithmic. The ±lcr regions around the rates from Okumura et al. 
(2014) have been traced by calculating the error using the statistical and systematic errors 
of the points measured by these authors. Note that the upper systematic errors include 50% 
extinction by dust, which explains why the upper error line is more distant from the main 
fit than the lower one. 







































22 


di^du) depends on the cosmological parameters Hq, Qm, and and 

P r^lim 1 

= TT / 7TT-i5^ + l)E]e(z,f!„,f!A)d2 (8) 

Hq Jo (1 + zy 

We adopt Hq = 70 km s“^ Mpc“\ VLm = 0.3, and = 0.7 in calculating our cosmic 
gamma-ray background. Indeed, the adopted values for VLm and Ua are the values favored 
by the PLANCK collaboration (Ade et al. 2013). The Hq = 70 km s“^ Mpc“^ value is a 
mean value coming from the discussions of the various approaches to determine Hq. 

Thus, the universe is flat and: 

e{z, LIm, = [(1 + zf(l + LImz) - z(2 + z)VLy\~^/‘^ (9) 

We take zum = 2.5, as our various test calculations show that the contribution to the 
background from SNe la at higher redshifts is negligible. That can readily be understood 
because SNe la rates drop signihcantly at redshifts 2: larger than ~ 2. Furthermore, in order 
to compare calculated fluxes with observed values, we need to divide the F^{E) above by 
dvr, to convert to the units of observed fluxes (photons cm“^ s“^ keV“^ sterad“^). 

In Figure 7, we hrst show the calculated cosmic gamma-ray background from SNe la 
model W7fm (W99b), for the SNe la rates of the Okumura et al. (2014). The rate is the 
main £t to the data of Okumura et al. (2014) at different redshifts 2; (the blue solid line 
in Figures 5,6). We show the predicted backgrounds when considering the contributions 
from SNe la up to 2; = 0.50, 1.00, 1.50, 2.00, 2.25, and 2.50, respectively (the increasing 
contributions correspond to the also increasing values of z). It can be seen that, indeed, 
after 2; = 2, further contributions to the background are negligible. In the Figure, the black 
squares correspond to the COMPTEL data, analyzed by Kappadath et al. (1996), while 
the continuous line gives the results from the Solar Maximum Mission (Watanabe et al. 
1999a), the dotted lines being the 1 a upper and lower limits. 


In order to make a consistency check, we calculate the gamma-ray background from 
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rates resulting from the fit to the D04 data (red dashed line in Figure 5 and Figure 6) and 
its Icr and 2 a predictions, using the W7fm model. Though the model is not the same as 
that used in AKH2005, we expect that our predictions should not signihcantly differ from 
theirs. Indeed, our resulting cosmic gamma-ray background, shown in Figure 8 as a black 
solid line, coincides with that of AKH05. In Figure 8, we also plot the la upper limit in 
blue dashed line and the 2cr upper limit of the cosmic gamma-ray background using the 
D04 SNe la rate as the red dashed line conhrming AKHOS’s conclusion that SNe la can not 
account for the observed gamma-ray background. However, that conclusion is mostly a 
direct consquence of the D04 SNe la rate used in the calculations, as we shall see below. 

Having calculated the gamma-ray background using the £t of the SNe la rate to the 
D04 data, we now consider the fit of Okumura et ah (2014). In Figure 9, we show the 
background spectrum for the mean value of the fit by Okumura et al. (the blue solid line in 
Figure 5 and 6) as the black continuous line. Clearly the best fit to the SNe la rate from the 
Subaru/XMM-Newton Deep Survey does not produce good agreement with the observed 
cosmic gamma-ray background from SMM and COMPTEL (Figure 9). In Figure 9, we 
also plot the background spectrum for the +- Icr limits of Okumura et al. (2014) shown 
in Figures 5 and 6 and their correspondings ±2cr bounds. We include the corresponding 
predictions in Fig. 9 for the gamma-ray background as they represent the most favorable 
case for SNe la as sources in the MeV range. 

In this work, for completeness, we calculate the background fluxes using the SNe la 
rates derived from the fit to the rates compiled by Graur et al. (2014) (whicch include 
their own CLASH data) (Fig 10) and by Rodney et al. (2014) (CANDELS) (Fig. 11). 
The Rodney et al. (2014) rates, shown as the green shot-dashed line in Figure 5 and 6, 
fall much lower than Okumura et al. (2014) rates. And the corresponding gamma-ray 
background fluxes are lower than the fluxes produced by using the rates from Okamura et 
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Fig. 7.— Contributions to the gamma-ray background when considering the mean value fit 
to the SNe la rates from Subaru-XMM (Okamura et ah (2014), for redshifts up to z = 0.50, 
1.00, 1.50, 2.00, 2.25 and 2.50 (curves in the order of increasing backgrounds at energy E 
around 100 keV correspond to the also increasing values of z), for model W7fm. One can 
see that, after z = 2, the additional contribution becomes very small. Here and in the next 
three Figures, the black squares correspond to the COMPTEL data, analyzed by Kappadath 
et al. (1996), while the continuous line gives the results from the Solar Maximum Mission 
(Watanabe et al. 1999a), the dotted lines being the 1 a upper and lower limits. 
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Fig. 8.— This figure shows the gamma-ray background obtained using the mean fit from 
D04 to the SN la rates (see Figure 4) as the black solid lines and the la upper and lower 
bounds are shown as blue dashes lines. The 2cr upper and lower bounds are shown as the 
red dashed lines. With those rates, the SNe la can not account for the observed gamma-ray 
background. We have used here the fully mixed W7 model (W7fm), as in all the other 
Figures. Indeed, our results for the cosmic gamma-ray background shown here as the black 
solid line coincides with that of AKH05. 
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Fig. 9.— This figure shows the gamma-ray background using the Okumura et ah (2014) 
fit to the observed SNe la rates (black continuum line), together with its ±lcr bounds (blue 
dashed lines) and ±2cr bounds (red lines). The ±la and ±2cr bounds have been traced by 
calculating the total error in the rates from the statistical and systematic errors of the points 
measured by these authors. The upper systematic errors include 50% extinction by dust. 
The model used is the fully mixed W7fm (W99b), but predictions for the models W7, W7dt, 
NlOO, and the two models of the violent merging of a couple of WDs are qualitative similar. 
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al. (2014) (see Figure 10, 11). HB2010 (Fig 12) prediction is shown here for the purpose 
of checking for consistency when reproducing their estimate. This estimate, based on more 
modest rates than those from Okumura et al. (2014) £t, would not give enough background 
to explain the MeV emission, either. 

One might argue that the SN la rates toward higher redshifts could be underestimated 
if one does not take into account the distribution of SN la light curve stretches, because 
that would disfavor the detection of small-stretch, lower luminosity supernovae. However, 
all the groups mentioned above have taken due account of this effect in modeling observed 
rates. 


6. Discussion 

SNe la are known to be prominent contributors to the annihilation line of positrons at 
511 keV and to the positronium continuum observed by OSSE an by INTEGRAL. Milne 
et et al. (2001) estimate that 30-50 % of Galactic positrons may be explained by SNe la 
and massive stars. The maximum of the contribution occurs at low bulge to disk ratios. 
Uncertainties of this contribution come from the SNe la Galactic rates (Milne et al. 2001). 
Such non-fully attribution to SNe la of the 511 keV line has led to the suggestion that 
most of the contribution comes from annihilation of dark matter (Fayet 2004; Boehm et 
al. 2004). In the GeV regime, the lack of a full contribution from astrophysical sources 
such as blazars, led to the proposal that weakly interacting massive particles (WIMPs) 
annhilations in the Galactic halo could be the main contributor to the diffuse GeV 
background fluxes. Those WIMPS would have masses in the range from 0.1 GeV to 10 
GeV, for a variety of dark-matter halo models (Pullen, Ghary, & Kamionkowski 2007). 
However, recent work (Ajello et al. 2015; Di Mauro & Donato 2015) greatly constrain the 
dark matter contributions in this energy range. These authors find that blazars, together 
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Fig. 10.— This figure shows the gamma-ray background using the Graur et ah (2014) 
compilation of SNe la rates (see Figures 5 and 6) as the solid line, and the upper la 
(blue dashed line) and 2a bounds (red dotted line). With those rates, the SNe la can 
not account for the observed gamma-ray background at the 2cr level (nor even at the 3 a 
level). The model used is the fully mixed W7 (W99b), but predictions for the other models 
are qualitatively similar. 








29 



Fig. 11.— This figure shows the gamma-ray background using the Rodney et ah (2014) 
(CANDELS) fit to the observational SNe la rates with their prediction (solid black line) and 
the upper 1 a (blue dashed line) and 2 a bounds (red dotted line). With those rates, the SNe 
la can not account for the observed gamma-ray background at the 2 a level (nor even at 
the 3a level). The model used is the fully mixed W7 (W99b), but predictions for the other 
models are qualitatively similar. 
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Fig. 12.— This figure shows the gamma-ray background using the HB2010 rates (private 
communication), together with their prediction (solid red line). The model used is the fully 
mixed W7 (W99b), but predictions for the other models are qualitatively similar. We use 
in red the calculations provided by HB2010. Our calculations using their rates are shown in 
black and they agree with theirs. 
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with star-forming and radio galaxies, can account for the gamma-ray background there. 

In a similar way, the early estimates of SN la rates at high 2 ; (D04), predicted a 
small contribution from Type la SNe to the gamma-ray background (AKH04). In that 
calculation, the SN la contribution is so small that even the 2 a upper bound derived from 
the D04 SNe la rates cannot account for the measured gamma-ray background (AKH04, 
and Figure 7 of this work). This opened a window for contributions from WIMP decays in 
the MeV range (see, for instance, AKH05; Cembranos, Feng & Strigari 2007), or from new 
populations of quasars (Ajello et al. 2009), amongst others. Here we show that the recent 
estimates of SNe la rates at high 2 , taken at face value, can not provide an explanation for 
the cosmic background in the MeV range. 

Recently, Ajello et al. (2009) analyzing the three year data from the Swift/Burst Alert 
Telescope survey concluded that the Flat-Spectrum Radio Quasars (FSRQs) could account 
for most, or even all, of the cosmic gamma—ray background for energies above 0.5 MeV 
by assuming its inverse Compton peak is located in the MeV band. We show in Figure 
13 as the long-dashed black line the best estimate of the contribution of FSRQs to the 
cosmic gamma-ray background. In Figure 13, the solid black line is the total contribution 
of FSRQs and SNe la when we adopt the 1 sigma upper bound of SNe la rates from 
Subaru/XMM-Newton (Okumura et al. 2014). The decline of the FSRQs which match the 
obeservations above 10 MeV is dN/dE ~ The shape of the SNe la contribution is 

very similar up to 3 MeV. Clearly the gamma-ray background from SNe la can improve the 
£t of ESRQs to the measured gamma-ray background between 0.3 and 0.8 MeV. 
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Fig. 13.— The cosmic gamma-ray background in the 1 to 30,000 keV range, from different 
sources. Shown in the dash-dotted dark green line is the predicted contribution by SNe 
la, when adopting the 2 a upper limit of the SN la rates from CANDELS (Rodney et ah 
2014), while in blue is that corresponding to the Icr upper bound of the rates by (Okumura 
et al. 2014) and in magenta we show the gamma-ray background from the upper 2 ct to 
the rates from Subaru/XMM-Newton (Okumura et al. 2014), the magenta line {2a upper 
bound) being very unrealistic. The long-dashed black line is the contribution from FSRQs 
(Ajello et al. 2009), while the solid black line is the total gamma-ray background of the 
long-dashed black line (FSRQs) and background from the for la upper bound of SNe I rates 
from Subaru/XMM-Newton (Okumura et al. 2014). This gives a good £t to the observed 
gamma-ray background of SMM (Watanabe et al. 1997) but only because of the FSRQ 
contribution. 







7. Conclusions 


Our calculations demonstrated that the use of SNe la rates ”as measured” at various 
redshifts are too low to account for the bulk of the observed gamma-ray background in the 
MeV range. We hnd that, above = 2-2.5, the gamma-ray emission from SNe la makes 
no signihcant contribution to the cosmic background. It is thus not expected that future 
determinations of the SNe la rates above that range would make a noticeable difference. 
However, what does make a difference is the overall SNe la rate from z = 0 to z = 2. We 
have shown here that recently reported values of the SNe la rates from Okumura et ah 
(2014), Graur et al. (2014) and Rodney et al. (2014) predict a contribution that is about a 
factor of hve below the observed background. Contributions from ^ beyond 2 are, as noted, 
too distant to be signihcant, even in the most optimistic cases, where rates continue to grow 
(see Figure 5 and 6). 

A relevant aspect of the uncertainties in determining rates at high z (> 1.4) is that the 
detection efficiency decreases rapidly with redshift, as observed bands shift farther into the 
rest-frame UV spectrum, where SN la emission is scarce (see the spectra shown in Riess et 
al. 2007). This causes the uncertainties to become very large at z ~ 1.6 (Barbary et al. 
2012). This question and the one of addressing dust obscuration which requires a model 
for cosmic dust evolution, makes the measurements of the rates at high z very unreliable. 
Note, however, the good agreement between the plethora of measurements of SNe la at low 
z (Figure 6). 

The current measurements of rates at high 2 : concentrate on normal (cosmological) 
SNe la. They make up to about 60% of all SN la explosions (see Ruiz-Lapuente 2014, 
for a review). The average predictions for the gamma-ray background, which at face 
values give gamma-ray backgrounds that are too low, might slightly change if the bulk 
of the non-cosmological SNe la were included. It has been estimated that up to 30% of 
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SNe la could be of the Type lax (Foley et al. 2013), though a lower estimate is found in 
White et al. (2014). Those SNe la produce smaller amounts of ^®Ni, thus making an even 
smaller contribution to the background (although their smaller ejecta mass increases the 
gamma-ray escape fraction). Furthermore, there is also a 15% group of SN1991bg-like 
events, which are much less luminous than the cosmological SNe la. 

From the opposite perspective, a factor that could make the mean gamma-ray fluxes 
larger would be to adopt models brighter in gamma-rays than W7 or W7fm (both models 
giving about the same average contribution over 600 days after explosion); some SNe la 
models can give twice as many escaping gamma-rays. One should acknowledge, however, 
that those models might fare worse in the comparison with observed spectra and light 
curves in the optical range. The hrst SN la observed in gamma-rays, SN 2014J (Churazov 
et al. 2014; Diehl et al. 2014), can be well reproduced within the set of realistic models 
presented in this paper. SN 2014J synthesized an amount of ^®Ni that falls within the 
standard range (about 0.5-0.6 Mq). That supernova is best explained as the explosion of 
a white dwarf close to the Chandrasekhar mass (Diehl et al. 2014). Having considered a 
wide range of models, we hnd that this is not the decisive point in regards to the question 
we address in the paper, i.e., whether SNe la are major contributors to the gamma-ray 
background in the MeV range. 

In summary, following many measurements and estimates of the SNe la rate as a 
function of redshift, we now have a more reliable estimate of the SNe la induced gamma-ray 
background in the MeV range than it was possible in the last XXth century and beginning 
of the XXIst. As the yields of different SN la models do not vary greatly, the conclusions are 
nearly independent of the particular model taken to represent normal SNe la. Attributing 
the gamma-ray background to SNe la appears tempting, as we know that they are 
significant y-ray sources, but within the current uncertainties SNe la do not seem able 
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to account for all of the MeV background. One needs to look for another source of the 
gamma-ray background in the MeV range, with FSQRs as a promising object class. SNe 
la can produce up to 30-50% of the gamma-ray background emission (if we take observed 
event rates at face value) and the remaining half to two thirds would have to be attributed 
to one or more contributors. 
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